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Abstract 

We discuss predictions for the lepton-flavour changing decays Z — > e/i, /ir, er 
which may be searched for at the Giga-Z option of the Tesla Linear Collider 
project with Z resonance production rates as high as 10 9 . We try to be as 
model-independent as possible and consider both the Dirac and Majorana mass 
cases. The Standard Model, if being minimally extended by the inclusion of light 
neutrino masses with some mixings as observed in neutrino oscillation search 
experiments, predicts completely negligible rates. With a more general neutrino 
content quite interesting expectations may be derived. 
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1 Introduction 



With the Giga-Z option of the Tesla project one may expect the production of about 
10 9 Z bosons at resonance 0. This huge rate, about a factor 100 higher than rates 
at LEP 1, allows one to study a number of problems with unprecedented precision. 
Among them is the search for lepton-flavour changes in Z decays: 



e/i, fxr, er. 



XI) 



Non-zero rates are expected if neutrinos are massive and mix 0, [| |J. Often one 
considers the branching ratio for the production of the following sum of charged states: 



BR(Z -> If If) 



T(Z -> hh + hh) 



X2) 



First predictions for flavour- changing Z decays in the framework of the Standard Model 
|| |, |7| , using techniques developed in || , were given in ||, [L0|, [TTj, 0] . The best direct 
limits are obtained by searches at LEP 1 (95% c.l.) 
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BR(Z 
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< 1.7 x 10~ 6 [14] 

< 9.8 x 10" 6 W% 

< 1.2 x 10~ 5 
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A careful analysis shows, taking into account realistic conditions at future experiments, 
that the sensitivities for the branching ratios could be improved considerably at the 



Giga-Z ||17|| , namely down to: 



BR(Z - 




BR(Z - 


-> e T r ± 


BR(Z - 


4 11*7* 



1-9 



< 2 x 10" 

< / x 6.5 x 10" 

< / x 2.2 x 10" 



;i.6) 

;i.7) 
:i.s) 



with / = 0.2-r 1.0. 

These numbers may be confronted with expectations derived from the signals for 
v n — v t oscillations in atmospheric neutrino experiments [18 , |I~9" , [2~0| , pi] , [22] , |23| , 24 1. 



They are at the 90% c.l. compatible with the following parameter set J2"5" 



Am 

Is 

sin 2 (2^ 



(2 4- 8) x 10~ 3 eV 2 , 
0.8 4 1. 



(1.9) 

;i.io) 



There is also evidence for u e — v„ oscillations from solar neutrino experiments EE 



, [3(J, |3l| , being compatible with: 

Am! „ - 



10 



-10 



10" 5 eV 2 . 



Xll) 



From reactor searches, there are no hints of v e — v T oscillations 32, 33 . For more 



details see e.g. the review [25 and references therein. 



2 



The good news from the evidences for neutrino oscillations is that they suggest 
non- vanishing rates for reaction (|1.1|) . The bad news is, that these rates are, if derived 
with ( |L9|) -( |ITTlT ) in the minimally extended Standard Model (z/SM), extremely smallQ: 

BR(Z /^r*) ~ 1(T 54 , (1.12) 
BR(Z e T ^) ~ BR(Z -> e^) < 4 x 10" 60 . (1.13) 

This is derived in Section |2| and is also in accordance with older calculations || [TT], [J3| . 

How do these small numbers arise? In Born approximation, the lepton-flavour 
changing Z decay into two charged leptons is forbidden in the uSM due to the GIM 



mechanism [36]. However, it may take place if n types of neutrinos have masses rrii 
and mix with each other (with mixing matrix Vy), i.e. if symmetry eigenstates and 
mass eigenstates are different in the lepton sector. Then, the virtual exchange of these 
neutrinos produces an effective lepton-flavour changing vertex and the corresponding 
branching ratio has the following structure: 

BR(Z-£J±) = — f 6 2 ^ |V(Ml)| 2 ^10- 6 |V(Ml)| 2 . (1.14) 

The form factor V(Q 2 ) depends on the details of the interaction: 

n 

V(Q 2 ) = Y, V ^ V ^/ M w)- (115) 
The vertex function V was calculated in 1982/83 independently by three groups || [10], 



IT] , [L2j1 for sequential Dirac particles, and in a more general context later |37|, p8| . 

The function V depends quadratically on the mass both in the small neutrino mass 
limit [3T5, [| and in the large neutrino mass limit [[5], [Tl], |TT], [T2] |: 



V(Xi < 1) - V(0) w (2.56 - 2.30 x z) + 0(X 2 i \nX i ), (1.16) 
7(Ai>l)-V(0) « i A. 4 + 0(lnA,), (1.17) 



with 



A. = ^. (1.18) 



The constant terms are not shown in (|1.16|) , ( |1 . 1 7|) since they drop out in (|1.15| ) due to 



the unitarity of the mixing matrix, and the branching ratio becomes proportional to the 
fourth power of the neutrino massses. It is this behaviour which makes the expected 
rates so extremely small for the experimentally evidenced tiny neutrino masses. For 
values of m 8 of the order of Mw, the vertex V is of the order one, and could become 
large if the rrii would be much bigger than Mw- 



1 Our estimate is in clear distinction to Eqn. (6) of p4|, where from the data a limit was derived 
which corresponds to BR(Z -> ^ t t ± ) « 0(1(T 8 ^ 10"^ 
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The evidence of tiny neutrino masses may also be indicative for a mechanism which 
produces at the same time very large masses. Heavy neutrinos are expected by some 



GUTs HU and string-inspired models [42, |43], |44J , and are suggested by the seesaw 



mechanism [ |5, |46|, |7| . Therefore, the above observations motivate us to have a closer 
look at the prospects of observing lepton-flavour changes with the Tesla Linear Collider. 
For the Giga-Z option, it might not be sufficient to apply the well-known and simple 
approximations for large masses m.j, but also the medium- or even small-mass cases 
may be of experimental interest. 

To be concrete, we will explore the following scenarios: 

(i) The z/SM. We treat the known light neutrinos (V e , z/^, v T ) as massive Dirac par- 
ticles. Individual lepton numbers L e , L M , L T are not conserved any more. The 
lepton sector is then in exact analogy to the quark sector. As a by-product, the 
Z decay amplitude into two quarks of different flavours can be read off from our 
general expressions. 

(ii) The uSM sequentially extended with one heavy ordinary Dirac neutrino. This 
case implies the existence of a heavy charged lepton as wellf]. It is not a very 
favoured scenario but we consider it as a simple application of the expressions of 
case (i) for heavier neutrinos. Again, total lepton number L is conserved. 

(iii) The vSM extended with two heavy right-handed singlet Majorana neutrinos. Not 
only individual, but also total L is, in general, not conserved since the presence 
of Majorana mass terms involves mixing of neutrinos and their charge-conjugate 
partners (antineutrinos), with opposite fermion- number. For two equal and heavy 
masses this case reduces to the addition of one heavy singlet Dirac neutrino |38f . 



In this latter case L is recovered ME 



With our numerical estimates we will be as model-independent as possible and will 
assume no constraints on neutrino masses or mixings, except for the ones imposed by 
the unitarity of the leptonic and neutrino mixing matrices, by the present bounds 
on lepton universality, CKM unitarity, and the measured Z boson invisible width 
49. 50, [51], |5j|. 54], and by oscillation experiments (see |25j for a review). 



In the following sections, we will discuss the predictions of the three scenarios for 
the lepton-number changing Z decay. In Appendices, the generalization of Lagrangian 
and Feynman rules of the Standard Model to the general case of Dirac and Majorana 
masses is explained, experimental limits on neutrino mixings and masses are quoted, 
and the calculation of the vertex function is sketched. 
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Dl: 




D4: vAAAAX,' 



DS: 



Figure 1: Feynman diagrams for the lepton- flavour changing Z decay. In the case 
of virtual, ordinary Dirac neutrinos, the ZviVj vertices in Dl and D3 are diagonal 
and the analogous quark-flavour-changing process can be obtained by replacing Ik by 
down-quarks and V{ by up-quarks. 

2 Predictions from the z/SM 

The amplitude for the decay of a Z boson into two charged leptons with different 
flavour, li and 1%, is given in a self-explanatory notation by: 



M 



I6irc w 



V{Q 2 ) £^ 2 (p 2 )7 M (l - i$)u h {-p 



where 



aw 



a 

2 5 
S W 



(2-1) 



(2.2) 



and the form factor V depends on Q 2 = (jp2 — Pi) 2 and can be written as: 



V(Q 2 ) 
V(\i) 



X)V W V^ V(A,), (2.3) 
i=i 

[v w (i) + vww{i) + V4,(i) + v H (i) + v W(j) (i) + , (2.4) 



with Vy being the leptonic CKM mixing matrix. In general, there are besides the vec- 
tor and axial- vector couplings fy and /a in also contributions of the fs, fp, /m> /e 
types, but for the production of on-shell fermions (with their masses being neglected) 
they vanish here. Further, it is fv = f a = V(Q 2 ) due to the presence of W bosons 
coupling only to left-handed fermions. The contributions from the one-loop diagrams 



2 A fourth generation of quarks is also needed to keep the theory anomaly free. 
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of Figure |T] depend on A* and additionally on 



A 



Q 



M 2 W 



which on the Z boson mass shell becomes 

\z - 



Ml 1 



M 2 r 2 
1VI W c w 



1.286. 



In terms of the usual vector and axial-vector couplings, 



Vi 
CI; 



I^-2Q i s 2 w ^If{l-As 2 w \Q i \l 



1 3 > 



the individual contributions are: 



(2.5) 



(2.6) 



(2.7) 
(2.8) 



From vertex diagrams: 
Dl : v w (i) = 

D2 : v ww (i) 
D3 : v^(i) 



D4 : v H {i) 
D5 : vw</>(i) 



+ a t ) [A Q (Co + C n + C 12 + C 23 ) - 2C 24 + 1] 
-{vi - at) Xi C , 

2c 2 w (2/^) [A Q (C n + C 12 + C 23 ) - 6C 24 + l] , 
A 2 

-fa + di) y c 



A; 



— (Vi — di 



Aq C 2 3 — 2C 2 4 + 



il-2s 2 w ) (2^)A,C 24 , 
-2s 2 w (2li L ) Xi C ; 



(2.9) 
(2.10) 



(2.11) 

(2.12) 
(2.13) 



From self-energy corrections to the external fermion lines: 

1 



DS: v s (i 



{vi Ac 2 w ai) [(2 + X i )B l + 1] . 



(2.14) 



The one-loop tensor integrals Co, Co, Cjj, Cy, and 5i are defined in Appendix |C[ With 
our numerical results for the Dirac case, we rely on two calculations, an old one |], [10 



and also this new, completely independent one. In the latter, the numerical evaluation 
of the tensor integrals is performed with the help of the computer program package 



LoopTools [55, 56 



Numerical results are shown in Figure |2|. The quantity presented is related to a 
branching ratio definition often used in the literature: 



T{Z - ft ft) 



2 x T{Z -> m) 

with N c as colour factor. Using 

r(z -> viv{) 



M? 



(2.15) 



(2.16) 



lcr 



7T J 16S 



Vi, Q 2 = M 2 Z 
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Figure 2: Contribution of one neutrino generation to the ratio B z of Eqn. ( \2.15] ) in the 
small and large mass regions for virtual, ordinary Dirac neutrinos and the analogous 
quark case. We set V^jV*^ = 1 here; dotted lines correspond to Q 2 = 0. 



we get a useful relation to the branching ratio introduced in (1^) and ( 1.14|) : 



BR(Z ^ ft ft 



B z = 0.1333 B z . 



(2.17) 



Figure |^ shows the contribution from one neutrino generation to the branching ratio 
as a function of the neutrino mass without the influence of the mixing matrix elements. 
We choose two interesting mass regions, namely that corresponding to the findings of 
neutrino oscillation searches (Figure 0(a)) and also that with potential predictions in 
the reach of the Giga-Z (Figure @(b)). The latter figure nicely agrees with the earlier 
calculations []T0] , |TT[ fl2]| , and the small mass limit with 0. The dotted lines in Figure 
^|(b) correspond to the approximation \z = 0: the large mass limit is reproduced quite 
well, while the small mass limit differs from the correct result. This is discussed in 



Appendix p. 3 



We will now estimate the branching ratio under the assumption that there are three 
generations of light neutrino flavours with a unitary mixing matrix V as evidenced by 
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experiment. The general form of this matrix may be chosen to be 





(Vel 


v e2 


v e3 


v,„ = | 










Vv T i 


v T2 


v T3 



C12C13 



C13S12 

i<5 












e i/3 









— C23Si2e — C 12 Si 3 S23 Ci2C2 3 e — S12S13S23 
■S23 s 12e l<5 — Ci2C2 3 Si 3 — C^^e* 15 — C23.S12.S13 C13C23 

(2.18) 

Here, we have three mixing angles and one CP- violating phase as in the quark CKM 
case, plus two CP-violating phases a,/3 if neutrinos are Majorana particles (they are 
strictly neutral so that less phase factors may be 'eaten' by redefining complex fermion 
fields). Current data suggest the following form of this matrix: 



C12 
1 



S12 





1 



1 r 1 
~V2 Cl2 V2> 



(2.19) 



where we have assumed |57[ a = (3 = 5 = 0, extracted from 0, £33| the si 3 = 0, and 
further assumed S23 = 1/ V2 (corresponding to maximal mixing) and left the #12 free. 

With 



T z = T(Z 
we get after trivial calculations: 



M w 
M z 
> all) 



80.41 GeV, 
91.187 GeV, 
2.49 GeV, 



(2.20) 
(2.21) 
(2.22) 



BR{Z 



lfl 2 , 



c? w M z 



V w V; 2 i[F(Ai)-V^(0)] 



1927r 2 c 2 y r z 

+ v, l2 vyv(A 2 ) - v(o)] + v, l3 v; 23 [v(A3) 



V(0)\ 



with 



a 3 w M z 



1.127 x 10" 



192 7 r 2 c 2 H/ r z 

For small neutrino masses, we show in Section p.2| : 

V(X i )-V(0) = a 1 (X z ) \. 



(2.23) 



(2.24) 



(2.25) 



The resulting branchings are, without approximations yet, but using the information 
from three-generation unitarity: 



BR(Z -> If if) 



u 3 w M z 
192 7 r 2 c 2 y r z 



\ai\ 



V W V* x Ai2 - V M n 3 A 



Z 2 3 A 23| , 



(2.26) 



with 

Ay = |Aj — Xj\ (2.27) 



and the ai is given in ( p. 67 ): 



|ai| 2 = 11.832. (2.28) 



There are two different cases to be considered when using now the specific mixing 
matrix (|2.19| ): 

c 2 s 2 

BR{Z -> e^) ss BR(Z -> e^r*) S3 1.333 x 1(T 5 \ 2 12 , (2.29) 

and 

BR(Z /^r*) = 1.333 x HT 5 - |s 2 2 Ai 2 - A 23 | 2 . (2.30) 

From the mass estimate ( |1.9| ) we get as additional input from atmospheric neutrino 
studies: 



(2 4- 8) x l(T 3 eV 2 



A 23 » V " ' ~ r ' = (3 + 12) x KT 25 , (2.31) 



1 vi- 
and from solar neutrino searches, ( 1.11|) : 

A12 sa 1.5 x (10~ 32 4 1(T 27 ) . (2.32) 

It is easy now to see that the expected rates for the lepton-flavour changing Z decays 
are limited to: 

BR(Z -> e^) S3 BR(Z -> eV) < 3.75 x KT 60 , (2.33) 
BR(Z -> ^r*) S3 (3 4 48) x KT 55 . (2.34) 

In Q2.33|) , we assumed arbitrarily a maximal mixing s 12 = l/y/2. These rates are 
extremely small. In fact, we will neglect the effects of the light neutrino sector in 
the next sections, where we extend the z/SM to accommodate heavy neutrinos, taking 
massless the known ones. 



3 Predictions from the z/SM plus One Heavy Dirac 
Neutrino 

Here, the only effective difference to the case before is the existence of a fourth gener- 
ation with a sequential Dirac neutrino of mass m^. In this case, the branching ratio 
gets: 

BR(Z-JfJ±) = °Vf Z r l V ^VU| 2 \V(\ N )-V(0)\\ (3.1) 
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with V(\n) given in Q2.4| ). The numerical results depend crucially on the mixing 
between the light and heavy leptons. An optimistic assumption would be maximal 
mixing: 

|V,„V i;K r=(-L) 4 = I. (3.2) 

This is of course unrealistic. Stringent (though indirect) limits may be derived from 
the analysis of flavour-diagonal reactions as advocated e.g. in |5(| 51], |52] as well as 
from the lepton flavour-changing process /x — > cy [pR A short summary may be 
found in Appendix [B] There the matrix B is, for this particular case of heavy Dirac 
neutrinos, B = V. 

In order to be definite, we show in Figure |3|a, solid line, the predictions for BR(Z — > 
/x ± r =F ) and assume the upper limit of the mixings allowed from ( |JJ.9| ) ( JR. 1 1| ) and ( |B.12| ): 

(V^V^I 2 < 1.5 x 10" 4 . (3.3) 

For the other two lepton flavour- changing Z decay channels, the corresponding graphs 
scale simply in accordance with the ratios of the mixing matrix elements: 

iVe^V^I 2 < 1.9 xlO" 4 , (3.4) 
(VeTvV^I 2 < 1.2 xlO" 4 . (3.5) 

These coupling factors are to be compared to ( |3.2| ): we observe a suppression of the 
expected branching ratio (for a given mass ttin) by more than three orders of magni- 
tude. 

As mentioned in Appendix [B], these bounds on light-heavy mixings from flavour- 
diagonal processes are improved by flavour- changing processes involving the first two 
lepton families. In fact, from [p0|1 

BR(fi -> ej) < 1.2 x 10~ u , (3.6) 

one obtains, for heavy neutrinos, the following (nearly) mass-independent limit J58]. |59| : 

|V e ivV; tW | 2 < 1.4 xlO" 8 , (3.7) 

much more stringent than ( |3.5| ). 

The Giga-Z discovery range is indicated in Figure |3]a, using the maximum values of 
the mixings allowed by ( |3.3|) . The range will be limited to the large neutrino mass limit 
if we believe in the relevance of the above mixing bounds. Then, the approximations 
of ( p.43|) apply. Numerically, this means for Xz = 1.286: 

\V(X N ) - V(0)\ 2 = - \ 2 N + 1.44 X N In X N - 3.49 \ N + 2.07 In 2 \ N + 0(\n\ N ). 

(3.8) 

If neutrinos with a mass of several hundred GeV or more would exist, there is a good 
chance to observe some effect of them from the Z decays under study. Indeed one 
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cannot constrain heavy neutrino masses from these processes since they also depend 
on the mixings, and vice versa. The fact that in Figure ^|a the discovery reach of LEP 
or Giga-Z cuts the curves means: if no event is observed, mixings must be smaller 
than the ones employed (we took present upper bounds) for neutrino masses above 
the intersection point; or if some effect was observed then the plot would provide a 
lower bound for the heavy neutrino mass. In fact, it is roughly only the product 
m % VjiivV* jy which can be constrained, assuming only one heavy Dirac neutrino. 



4 Predictions from the z/SM plus Two Heavy 
Right- Handed Singlet Majorana Neutrinos 

Some basic features of Lagrangians with Majorana mass terms and their relations to 
the simpler Dirac case are summarized in Appendix [A|. For a derivation of Feynman 
rules with Majorana particles |HJ we refer to The couplings of the virtual 

neutrinos to the W bosons are left-handed, v i + aj = 1, Vi — cii = 0, and 2Jg £ = 1, 
while those to the Z boson and to the Higgs particles are non-diagonal and contain 
right-handed admixtures. This may be seen from ( |A.44| )-( |A.46| ). The amplitude for 
the decay (|1 . 1| ) is again given by Eqn. ( |2.1|) , but this time the form factor Vm is non- 
diagonal not only in the external charged leptons, but also in the virtual neutrinos due 
to the ZviVj coupling, see ( |A.45|) : 



V M (Q 2 ) 



n G +n R 

J2 B W B* . V(i,j), 

V(i,j) = V(\i, Xj, Cij) 

= [vw(i,j) + 8ijV W w(i) + v^ij) + SijV^i) + SijV W cf>(i) + <%^e(0] • 



(4.1) 



The new non-diagonal terms arise in diagrams Dl and D3 of Figure 0: 

Dl:v w (iJ) = -C i:j [X Q (Co + C n + C X2 + C 23 ) - 2C 24 + 1] 
+C*j a/ AjAj C , 



D3 : Vj(i,j) 



~TT~ + V 9 — 



Xq C23 — 2C 



2d 



(4.2) 

(4.3) 
(4.4) 



The matrices B and C are introduced in Appendix [A.3| . Both contributions are quite 
similar to ( |2.9p and ( |2.11|) when there the right-handed couplings are retained. The 
other contributions are the same as in the ordinary Dirac case. The vertex reads: 



n G +n R 



V M (Q 2 



3,^. [^•F(A i ) + C ii C(Ai,A i ) + C*.V^iy(Ai,A^)J,(4.5) 



to be compared to the case of nn heavy sequential Dirac neutrinos (B = V) : 

nG+n H 

V(Q 2 ) = BiiiB* 2 j V(Xi). (4.6) 



i=l 
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Figure 3: Predictions for the branching fraction BK(Z — > \xt + fif) as defined in 
the text. Indirect upper limits on the mixing angles from flavour- diagonal observables, 
obtained in the models used, are taken into account. 



We now have to go into a specific model in order to make definite predictions. For 
the case chosen, namely the Standard Model extended by two heavy Majorana singlets, 
the sums over virtual neutrinos involve nc=3 light and n^=2 heavy neutrinos. The 
model is described in Appendix [A.3| . For some important properties of the mixing 



matrices B and C see ( |A.49| )-( |A.52j ). The light neutrino sector is practically massless. 
Using the properties (|A.49Q - (|A.52|) , one can then write the form factor Vm in terms of 
the heavy-neutrino sector only. Actually, it is: 



V{X i ) = F{X i ) + G(X i ,X l ). 



(4.7) 



Further, the functions F, G, H are uniquely defined and from (|A.49|) -( |A.52| ), and taking 



A, = for i 



BR(Z -> If l d 2 



,Ug, it is straightforward to prove that: 



c? w M z 



192n 2 c 2 w r z 



x {5 NiNj [F(X Ni ) - F(0) + G(X Ni , 0) + G(0, X Ni ) - 2G(0, 0)] 
+C NiNj [G(X Ni1 X N .)-G(X Ni ,0)-G(0,X N .) + G(Q,0)} 



+ C *N t N, \ ^n,X Nj H(X Ni ,X Nj )} 
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It is possible to express the couplings B, C on the right-hand side by the mass 
ratio r = m 2 N2 /m 2 Ni plus the three light-heavy mixings s Ve , s v , s Vt . The relations 
are explicitly given in ([A.53|) — ( |A.57 ). The upper limits for these mixings (given in 



Appendix |B.1| ) have been used to obtain the graphs in Figure 0(b). The mass ratio 
r has been taken as a free parameter. Perturbative unitarity constrains the masses of 
the neutrinos so that they cannot be arbitrarily heavy. This is discussed in Appendix 



B.2| . We see again from the figure that Giga-Z has a discovery potential, preferentially 
in the large neutrino mass region. Similar curves can be obtained for Z — > e/i and 
Z — > er. 

The distinguished case of two Majorana singlet neutrinos with equal masses, forming 
effectively a singlet Dirac particle ||38|| , results for r = 1, see the solid line. This line 
has been taken over from Figure 0(a), there as the dashed line, in order to show that 
due to the different coupling structure the simple sequential Dirac neutrino case does 
not constitute a limiting case for large masses. The deviations are due to the terms 



from the non-diagonal elementary Z couplings proportional to C and C* in (O). The 
C* terms drop out for r = 1. 

In contrast, predictions for Majorana and Dirac neutrinos approach each other in 
the limit of small masses. This may be seen from (18) using the unitarity relations 



( |A.49| )-( [A.52| ) and the Taylor series expansion of the vertex function in powers of the 



neutrino mass (Appendix [D . 2|) . This phenomenon is just another example of what is 
called in the literature the "practical Dirac-Majorana confusion theorem" |6l| (see also 
the recent discussion in |6C| |67j and references therein). 

At the end of this section, we would like to comment on the literature for the 
reaction ( |1 . 1| ) . The early papers did not include the more "realistic" and physically 
most interesting models with Majorana masses. These cases were studied in detail by 
38| and, in the context of left-right symmetric models, by fB5 |. While we reproduce 



the large mass limit given there, we obtain slight deviations in the medium mass case 



as given in Eqn. (Bl) of [3S]. Further, the Figures 7 and 8 in [38] are not exactly 
reproducible in the intermediate mass range. Bowever, the large mass limit seems to 
be the only potentially relevant case. 



5 Summary 

From our study of the decays Z — > e/i, er, /it, in the context of the Giga-Z option of 
the Tesla linear collider project, we conclude: 

• Neglecting the influence of the mixing angles, the expected branching ratios de- 
pend on the neutrino masses m ; and are of the order BR ~ 1CT 5 (mj/Mw) 4 both 
in the small and large mass limits. 

• If there exist only the known three generations of ultra-light neutrinos, m v <C 
M w , there is absolutely no hope to see any effect. 
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• For heavier neutrinos with masses of the order of the weak scale, m u ~ Myy, one 
has to calculate the form factors describing the vertex without any approximation. 
The necessary exact formulas have also been given. 

• The light-heavy mixing angles are not so strongly restricted as with the minimal 
(one-family) seesaw mechanism when interfamily seesaw type models are consid- 
ered. But, unfortunately, these mixings have already been constrained to be very 
small, so that the lepton flavour- violating reactions under study are beyond the 
discovery reach of the Giga-Z for m v ~ My/ . 

• Since, ignoring the light-heavy mixings, the branching ratios are proportional to 
the fourth power of the neutrino masses, there is a discovery potential in the large 
mass case, m v ^> Mw- Nevertheless, they are limited in practice by potentially 
small mixing factors (bound independently of the heavy neutrino masses by other 
experiments) and upper neutrino mass bounds from unitarity considerations. 

• In fact, there is an interplay between heavy masses and light-heavy mixings: 
the mixings must be small for very large neutrino masses, since otherwise the 
scattering matrix elements would grow above the unitarity limit. 

Summarizing, the Giga-Z offers nearly three orders of magnitude gain of sensitivity 
compared to LEP 1. This opens quite interesting opportunities to search for lepton- 
number changing processes, if there exist heavy neutrinos sufficiently mixing with the 
light sector, within a quite broad allowed region according to present limits. 
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A Lagrangians and Feynman Rules 

We make extensive use of the notion of Majorana particles, since in general neutrinos 
may be of this type, and in fact in GUTs often exactly this happens. Majorana particles 



fl6l| are neutral fermions ip, fulfilling 

r = *p, (a.i) 

where ip c is the charge-conjugate of i/j. 

Some introduction on notations are given in Appendix |A.1| , We observe experimen- 
tal evidences for ultra-light neutrinos in the known fermion families on the one hand, 
and on the other there are unifying theories with potentially ultra-heavy Majorana 
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neutrinos. That both phenomena might be related will be made plausible by a toy 
example in Appendix |A.2| . A more general ansatz for Majorana mass terms is finally 
given in Appendix |A.3 . 



(A.2) 



A.l Dirac neutrinos rewritten 

The mass-term Lagrangian for Dirac neutrinos is: 

- C D = m D {~ L v R + - R v L ) ee i (F) "J") , 

where we introduce (self-conjugate) Majorana fields (%°) 

m - ft: 3). 

with 

z/ c ee Cz/ T (A.4) 

and C being the charge-conjugation matrix. The mass matrix can be brought into a 
diagonal form in the basis of (x): 

- \ @(7 m ° D ) «) = i". D (66 + 66). (A.5) 

In the last step, the field component £1 is introduced as a chiral transform of Xi i n 
order to make the positive mass eigenvalue explicit. The fields (x) and (x°) are related 
by the unitary matrix U, and £1 additionally by 75: 



(A.6) 



6 = 75Xi = -^(-vl + vI-vr + Vr), (A.7) 
6 = X2 = -^=K + z/£ + ^ + ^). (A.8) 

With the above chain of equations we have shown that one ordinary Dirac neutrino 
is equivalent to two Majorana neutrinos of equal mass and opposite CP parities. Ev- 
idently, if the mass matrix is not of type ( |A.2|) , then true Majorana particles are 
realized. 
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A. 2 The seesaw mechanism 



The seesaw mechanism [45, 4J, 47] allows to understand the lightness of the known 
neutrinos by the introduction of heavier ones. However, it is not the only possible 



solution to this puzzle; see e.g. |59], |69| for a nice discussion. 

In Appendix KA. it was shown that a Dirac-neutrino mass term may be written as 
symmetric 2x2 matrix with only off-diagonal elements. Consider now the case where 
only one right-handed singlet neutrino is added to the ordinary left-handed doublets 
of the SM. Let's take one generation of left-handed light neutrinos for simplicity and 
with no loss of generality. The corresponding mass matrix is in general 



M 



and can be diagonalized by 



U 



m L m D 
m D m R 



cos 9 V sin 9 V 
- sin 9 V cos 9 V 



(A.9) 



with 



tan 26^ 



Ivtid 



cos 29 



m R - m L 



a/ ( m R ~ wil) 2 + 4m 



m R - m L 

yielding two eigenstates with different masses 

^ jmjj + m L =f [{m R - m L f + Am 2 D ] 1/2 J 



(A.IO) 



(All) 



m u ,m N 



(A.12) 



That is, the most general mass term of one four-component self-conjugate field describes 
two Majorana particles with different masses. There is a particular configuration that 
corresponds to one Dirac neutrino, as shown in Appendix A.l . 

In the SM, m L = since the SM Higgs sector consists of a Higgs doublet (see 
Appendix |A.3| ). Take now m R ^> mp. Then, the two physical states^ are a light and 
a heavy neutrino with masses 



m N 



ml/m R , 
m R , 



and the light-heavy mixing angle is: 

s u = sin9 u ~ m D /m R ~ \/m v jm N . 



(A.13) 
(A.14) 



(A.15) 



This is the minimal seesaw mechanism: the larger the heavy scale is, the smaller the 
light neutrino masses are, and, at the same time, the smaller the mixing becomes. 



3 Actually, one of the mass eigenvalues in (A.12) is negative. This is not a problem since then the 
true mass eigenstate is a chiral transform of the original field which has a mass term of opposite sign 
as in (TO) fo|. 
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Actually, by taking mjy > Mz and the present bounds on the light neutrino masses 
0, 





< 


2.5 eV, 




< 


160 keV, 




< 


15 MeV, 



(A.16) 
(A.17) 
(A.18) 



we get from (A..11) the light-heavy mixing angles: 

-ii 



< 


< 


3 x 10 


< 


< 


2 x 10 


< 


< 


2 x 10 



(A.19) 
(A.20) 
(A.21) 

These mixing angles are too small to be constrained by the experimental LEP and 
low-energy limits @ 0. 

Nevertheless, one might have ^ 0, by introducing a Higgs triplet in the SM, [] 
and hope for a conspiracy between and to get light-heavy mixings of order 
one. This does not seem very reasonable. Thus, the minimal seesaw apparently lacks 
of phenomenological interest. 



A. 3 Majorana neutrinos 

In the interfamily seesaw-type models, several right-handed neutrinos are introduced 
and large light-heavy mixings can be obtained even for = |7T|, [72|, |73| . 

Therefore following [38], we now introduce in our study two heavy right-handed 



singlet neutrinos and treat the light-heavy mixings s Ul and the heavy masses , m^ 2 
as free phenomenological parameters. 

Consider an extension of the SM that incorporates right-handed neutrinos z/r 
(singlets under SU(2)(g>U(l)) to the already present set of n G = 3 generations of left- 
handed neutrino doublets. We keep the Higgs sector untouched. 

It is convenient to arrange all the independent neutrino degrees of freedom into two 
vectors of left-handed fields: 

(z/°) = (u e ,u M u T ), (A.22) 
(N°) = 04), z = l,...,n fl . (A.23) 

The charge-conjugates of right-handed neutrinos have been introduced: 

v% ee Cv£. (A.24) 

The conjugate field has chirality and lepton/fermion numbers opposite to the original 
field (e.g. the z/^'s are left-handed). 



4 The introduction of a Higgs triplet affects the value of the parameter p = / \M\c 



w ) 
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In the basis 

(n°) = (y°, N°) (A.25) 
the most general mass-term Lagrangian reads: 

-C M = l(ri^)M(n )+h.c. = hn ) T CM(n°)+h.c., (A.26) 
where the Majorana fields 

( X °) = (n°) + (n oc ), (A.27) 

containing both chiralities, have been introduced. 

The mass matrix M is symmetric and can be written in a block form: 

M= /m L m£\ 
^ m D m R J 

It has dimension + Ur and can be diagonalized by a unitary matrix U of the same 
dimension, 



M = U T MU = (U*) f MU = dia.g{m Ve ,m Vll ,m UT ,m Nl ,m N2 ) 

~ diag(0, 0,0,m Nl ,m N A, (A.29) 

since we are mostly interested in a heavy neutrino sector consisting of two right-handed 
neutrinos. The interaction eigenstates, in terms of the left- and right-handed compo- 
nents of the mass eigenstates, are given by: 

Left-handed: (n°) = f „ ) = f Vl c ) = U P L ( X ) (A.30) 



,o ( - 



Right-handed: (n° c ) = ( ^ QC ) = ( ^ j U P H (\). (AM.) 



with 



P r ,l = \{±±1s). (A.32) 

The diagonal blocks in M, connecting states with opposite fermion number (m^, vcir) 
are called "Majorana" mass terms, while the off-diagonal ones (nip) are "Dirac" mass 
terms that conserve fermion number.^ The mass terms arise, after the spontaneous 
symmetry breaking (SSB) of SU(2) L (g)U(l)y — >U(1)q, from the Yukawa coupling of 
the fermion fields to the neutral Higgs field. The SU(2)l®U(1)y quantum numbers of 
a Higgs doublet and the fermion bilinears that can be constructed from the doublet 
and singlet neutrinos are 

*=($o) ~ (2,1) (A.33) 

AL= 0: V^v L ~ (1,0)® (2,-1) = (2,-1) (A.34) 

AL = ±2: 17[v L ~ (2,-1) <g> (2,-1) = (1,-2) © (3,-2) (A.35) 

AL = ±2: ~W R v R ~ (1,0) <g> (1,0) = (1,0) (neutral). (A.36) 

Non-diagonal Dirac mass terms produce transitions between states with different individual lepton 
numbers but the total lepton number, i.e. the fermion number, is conserved. 
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Therefore in a model with only Higgs doublets, the entry = 0. 

Our convention for the covariant derivative acting on a fermion doublet is 

D M = d, + ig T - ■ W, + ig'^Bp (A.37) 

with 

Q f = l{ L + Y/2. (A.38) 



After SSB one gets 



e = gs w = g'c w , (A.39) 

(A.40) 



Z»\ = f c w -s w \ ( Wl 
A? J \ s w c w J \ Bf, 

= ±={Wl±iWl). (A.41) 

The interaction Lagrangians of neutrinos with W, Z and charged Goldstone bosons 
(j) , in the weak basisf] 



C 



n G 

9 



ir 



-^W-^^P^+h.c, (A.42) 

i,.7=l 

9 
2cw 



"G 

= ~^Z,J2 frT k rPLVL k -W k rP R v C L k ], (A.43) 
fe=l 

become, in terms of physical neutrinos Xj with masses rrij, j = 1, no + tir: 

ug nG+n R 

C w = E BukfPLXi + te., (A.44) 

»=i j'=i 



£z = -ir-Zv E ^(^Cy-P^OXi, (A.45) 

n G n G +ra H , v 

^ ^ i'-E E ^ ~ ^) X, + (A.46) 

2 1 J 1 



with 



>l G 



B ki = E r w% (A.47) 

fc=i 

"G 

fc=l 

6 The interaction of the Z to two (right-handed) charge-conjugate neutrinos (non-existent in the 
ordinary SM) has been included in ( |A.43| ), making use of the bilinear transformation properties under 
charge-conjugation: ^Pj^Pl^iP — * —^I^Pr,^ - 
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V^j and By are the leptonic CKM mixing matrices and its generalized version, re- 
spectively |62|, 7^]. For Dirac particles B = V, Cy = 5ij and C*j = 0. For Majorana 
particles, in contrast, there are NC couplings of different flavours (FCNC) and with 
both left- and right-handed components. 

The matrix Vy is quadratic of dimension nc and By is rectangular no x (no + ^ij) 
and incorporates lepton-flavour changing mixings. The matrix Cy is quadratic, has 
dimension (nc + %) , and causes flavour non-diagonal Z\iXj interactions. 

These interaction Lagrangians involve Majorana fermions that complicate the eval- 
uation of S matrix elements, since extra Wick contractions survive in comparison to 



the case with only Dirac fermions. Following [74| [64| , one can write Feynman rules re- 
sembling the ones of the Dirac fermions, based on the well defined fermion flow, rather 
than on the fermion number flow which is not preserved in the vertices with Majorana 
fermions: after fixing an arbitrary orientation (fermion flow) for a given diagram, the 
vertices can be read off from the Lagrangian as usual, but for every vertex f\Tf 2 one 
has to add the reversed one, /fr"/|, with V = CY T C~ l . The fermion propagators are 
the usual ones. This effectively yields the same result for the vertex WUxj but a factor 
two larger for the vertex ZxiXj i n comparison to the case of Dirac neutrinos, since for 
two Majorana fermions (x = X c ) T = r". 

Important note: elsewhere in the text we refer to and N± as the neutrino physical 
states, rather than Xh to simplify the presentation, but with no possible confusion, 
since we always work in the physical basis. 



The matrices B and C obey a number of useful relations |37j : 



riG+riR 



B *ii B *2j 

/j C ik C jk 



fillh, 
n G 



(A.49) 
(A.50) 



fc=i 



k=l 



na+n R 



/ J B^Cfcj = B 



k=l 
n G +n R 



riG+n R 



k=l 



k=l 



k=l 



(A.51) 



} u m k C ik C jk = ^2 m k Bi k C* ki = m k B hk B hk = 0. (A. 52) 



In the case of hr = 2 the matrix elements involving heavy neutrinos can be obtained 
from ( A.50 ), ( A. 52] ) in terms of the light-heavy mixing angles (BJ3) and the ratio of 
the two heavy masses squared r = m% 2 /m 2 Ni , assuming that the light sector consists 
of massless neutrinos f38| : 



,1/4 



B 



l k N 2 



Vl + r 1 / 2 
i 

Vl + rV2 



(A.53) 
(A.54) 
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and 



J N 2 N2 



r l/2 "G 

l-7I7l£<> (A.55) 

fc=i 

lT^E<' (A.56) 

fc=i 

a/4 ™g 



WxN 2 ~ ~Cn 2Ni - , 1/2 J] S l ( A - 57 ) 

fc=l 



B Constraints on Heavy Neutrinos 
B.l Experimental bounds 

Several kinds of constraints on heavy-neutrino masses and light-heavy mixings can be 
obtained from experiment. 

Direct production searches establish the following limits on the neutrino masses at 
95% c.l. @: 

Stable neutrinos: m N > 45.0 (39.5) GeV (B.l) 
Unstable neutrinos: m N > 69.0 (58.2) GeV, (B.2) 

for Dirac (Majorana) particles, respectively. 

A general formalism to describe light-heavy mixings was developed in |49|, [58] . The 
mixing angles in our notation |5S| correspond to 

£l B ^.| 2 ' (B.3) 



Indirect constraints on the masses and bounds on the mixings are provided by two 
categories of LEP and low energy experiments: 

(i) Flavour-diagonal processes. They include mass-independent and model-indepen- 



dent light-heavy mixing constraints [49. 50, 51. 52] from tests of lepton universality and 



CKM unitarity and measurements of the Z boson invisible width, as well as other less 
sensitive studies like W mass measurements and low energy experiments like neutrino 
scattering, atomic parity violation, etc. Flavour-conserving leptonic decays Z — ► 
depend on masses and mixings through loop contributions and provide alternative 



constraints 53 



^ii) Flavour-changing processes. They include rare processes like \i —>■ ej, \x 



ee + e~, ji — >• e conversion in nuclei, r — > l a l^l~ and Z — ► ||9], |59|. They are 

mass- dependent (except for \i — » ej to a good approximation). One can get also less 
stringent light-heavy mixing constraints from oscillation experiments |58|, |5I| . 

The most stringent present bounds on the light-heavy mixings are provided by 
the flavour-diagonal processes. Exceptions are the ones involving the first two lepton 
families such as fi — » cy, Z — > e/i and \i — > ee + e~ 
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For illustration we show an example of the flavour-diagonal constraints. The effec- 
tive muon decay constant is related to the coupling Gf of the standard model by 



G, 



c„c,,,,G p. 



(B.4) 



The unitarity constraint for the first row of the CKM quark-mixing matrix implies |5T] 



Cv e Gl 

G,. 



— = 0.9992 ±0.0014. 

cr, 



:b.s) 



As a summary, from [[49], [5(], |5T|] (90% c.L): 



sl e < 0.0071 (0.005), 
< 0.0014, 



si T < 0.033 (0.01), 



(B.6) 
(B.7) 
(B.8) 



where the most conservative bounds are obtained assuming any kind of heavy neutrinos 
and the ones in brackets correspond to the case of SU(2) singlets. 

From the most recent update[| by |52|, assuming only heavy singlets, one gets (90% 
c.L): 

(B.9) 
(B.10) 

(B.H) 



< 


< 


0.012, 


< 


< 


0.0096 


< 


< 


0.016. 



A final remark is in order here. Using the Schwartz inequalities [pS 

iE B 



2 <s 2 u s 2 vl 



(B.12) 



one can infer indirect upper limits on the off-diagonal mixings (relevant for the flavour- 
changing processes) from the previous flavour-diagonal constraints. Nevertheless, for 
our scenario (iii), as already mentioned, the mixings can be obtained exactly from the 
properties of B and C and such inequalities are not needed. 



B.2 Decoupling and neutrino-mass upper limits from pertur- 
bative unitarity 



The heavy-neutrino masses are restricted by the perturbative unitarity condition on 



the decay width of heavy neutrinos [75, 76, 77 



Ni 



(B.13) 



7 These latest bounds are more conservative than the earlier ones in |5lJ] due to the fact that 
present determinations of the elements of the first row of the CKM matrix are not compatible with 
unitarity, and hence this constraint is eliminated from the analysis. 
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The total decay width of a heavy Dirac neutrino (four d.o.f.) with mass 3> 
M W ,M Z ,M H is iJ, 0, p| : 

3 3 



AT* 



J]r(iVi -> z- w+) + ^ [T{Ni -> ^ fe z) + r(iv 4 -> ^#)] 



fc=i 



/t=i 



|2 

l k Ni\ , 



(B.14) 



fe=i 



and a factor two larger for a heavy Majorana neutrino (two d.o.f.). Therefore, the 
perturbative unitarity bound expressed in QB.13|) reads 



m 2 Ni 



Em. 



m 2 AT . C 



< 



k=l 



4M w /aw for a Dirac neutrino 
2M w /aw for a Majorana neutrino, 



(B.15) 



(no summation over repeated indices is understood) which shows implicitly that heavy 
neutrinos decouple [|7|, |79fl , in accordance with the Appelquist-Carazzone theorem |8l| , 
p2| : the unacceptable large-mass behaviour of the amplitudes (~ m 2 N ) is actually cured 
when the light-heavy mixing (~ m~^ 2 ) is taken into account |^3|, Qj. 

Taking the values of in terms of the light-heavy mixing angles ( |B.3| ) for one 

Dirac or for two heavy Majorana neutrinos ( |A.53|) , (|A.54|) one can get the following 
upper limits: 



2 . 4M^ 

m N < 

aw 



-i -l 



E 

.k=l 



(B.16) 



for a heavy Dirac neutrino, and 

1 



2^-2 ^ 2M w l+r^ 
m N = -m N < 

r 2 a w r 



E 

,fc=i 



-i 



(B.17) 



for two heavy Majorana singlets. The latter bound is very stringent when and 
m^ 2 are very different. The upper mass limits on heavy neutrinos, from ( [B.16D , ( |B.17| ) 
are then: 



scenario (n): 



m 2 N 



< 



(4.2,4.4 TeV) 2 



scenario (m): m N ^ = z m N 2 



1 + r 1/2 

< — x (3.0,3.1 TeV) 2 , 



using the bounds (p.6j )-( fBT8|) and (|B . 9|j — ( |B ■ 1 1|) , respectively. 



(B.18) 
(B.19) 
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C The Vertex Function 



We use the notations of [Sl|. The vertex function contains the following one- loop 
integrals in D dimensions || 36 : 



16tt 2 



A{ml) 



A-D 



167T 2 



4-D 



16tt 2 



with 



and 



{Co,C^Cn(piQ,Pi;<rnim 2 2 



Q 2 = (P2-Pif 



A-D 



d D q 1 

d D q {l,g"} 
(2tt)^ £> £>i ' 

(2vr)^ XTOI> 2 



(C.l) 
(C.2) 
, (C3) 

(C.4) 



V 

v 2 



2 2 i • 

g - m + ze, 

(g + pi) 2 — m\ + ie, 



(Q + P2) 2 - "ml + 

They are decomposed into tensor integrals according to their Lorentz structure: 



p^B 1 , 



12 j 



ptPiC 2 i + p^c 22 + ip^vl + P%Pl)C 23 + g» u C 24 . 
We employ the following dimensionless tensor integrals and their abbreviations: 



(C.5) 
(C.6) 
(C.7) 



(C.8) 
(C.9) 
(CIO) 



Br 


= Si(0;Ai,l) 






= B^mlM^), 


(C.ll) 


C{0, 11, 12,21,22,23} 


= C{...}(0, Aq, 0; 1, Aj, Xj) 






= M^C { ... } (0,Q 2 ,0;M^,m J 2 ,m J 2 ), 


(C.12) 


C24 


= 6*24(0, Aq, 0; 1, Aj, Aj) 






= C 24 (0,Q 2 ,0;M^,m 2 ,m 2 ), 


(C.13) 


C{0,11, 12,21, 22,23} 


= C { ... } (0,A Q ,0;A i ,l,l) 






= Af^C { ... } (0,Q 2 ,0;m 2 ,M^,M^), 


(C.14) 


C*24 


= C M {0, \ Q , 0; A<, 1, 1) 






= C 24 (0,Q 2 ,0;m 2 ,M^,M^). 


(C.15) 



On the Z mass shell it is Xq = \z- The Aq, A^, Aj are introduced in 
(I2T 



and 
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The explicit expressions for the loop functions are: 



C{0,11,23} 

Ci2(Aj, Xj) 

C*24 



dx 



dy 
V 



1 



2(D-A) 2 



dx / dylnV, 



(C.16) 
(C.17) 
(C.18) 



with 



V 



Di 



ijk 



X Q xy + (X k - Xj)x + (-X Q + X { - X k )y + Xj - ie, 



(C.19) 



and the correspondences are: for the non-abelian diagrams (with elementary ZWW, 
ZW<p, Zcfxf) vertices), (i,j) are virtual W, <f> bosons and k a neutrino, while for the 
abelian diagrams k is the W,4> boson and (i,j) are neutrinos. On the Z boson mass 
shell: 



V 
V 



D ijW = X z xy + (1 - Xj)x + [-Xz + (Xi - l)}y + Aj - ie, 

ie. 



D W wi = X z xy - (1 - Xi)x + [-X z - (Xi - l)\y + 1 
In the Dirac case, it is A, = Xj. Further, 

1 



(C.20) 
(C.21) 



D-A 
1 



/ xdx \n[x + Xi(l — x) — ie] 
Jo 



X; 



D-A 2(1 -A,) 



1 + 



Xj In Xj 
T^A~ 



1 

4' 



(C.22) 



As may be seen, the tensor integrals Bi, C24 and 6*24 are ultraviolet-divergent in D 
dimensions. We mention that the functions are defined here with Minkowskian metric, 
so that all the functions introduced, except for B , Bi, C24, C24, have different sign from 
those used in e.g. [|TI| (Euclidean metric). In principle, the divergent parts depend on 
the regularization scheme and could differ by a finite, universal term yet. Recalling the 
UV-behaviour of the divergent integrals Bi, C24 and C*24 (see flC.22|) and ( |C. 18|) ) , we 
get divergent, mass- dependent contributions from individual diagrams to the vertex V. 

For the Dirac case: 

~1 A,: 



v^(i) 



(vi + a«) x [finite] 



a* x 



(2If) (I -2s 



XV) 



1 A, 



2D 



2D -A 
+ finite 



+ finite 



~ -(vi + di - Ac w ai) 



Xi 



D-A D-A 



+ finite 



(C.23) 
(C.24) 
(C.25) 



The sum of the terms proportional to Xi/(D — A) vanishes. The constant divergent 
terms (not shown here) sum up for individual vertices V(Xi), but vanish due to the 
unitarity of the mixing matrix for the complete vertex function ( |1.15| ): 

1 



V 



(2L 



D-A 



S, 



hh 



0. 



(C.26) 
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For the Majorana case {2P 3 L = 1, V{ = ai = 1/2): 



v<M>(i,j) 



-Cij x [finite] + C*- x 



5 in -( 



2s z w ) x 



1 yAfAj 

2 D-A 

A 



+ finite 



2D-4 
2 A, 



+ finite 



L>-4 D — A 



finite 



(C.27) 

(C.28) 
(C.29) 



The mass- dependent divergent terms in ( |C.28| ) and ( |C.29| ) cancel each other, and the 
one proportional to C* 3 - in (|C.27|) drops out due to the unitarity condition ( |A.52|) when 
the sum over virtual neutrinos is performed in ( j4.1|) . The constant divergent terms 
vanish again due to the unitarity relations ( |A.49| )- (|A.52| ): 



V 



M 
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\v — 4 ll12 



(C.30) 



D The Vertex Function for Large and Small Neu- 
trino Masses 
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D.l The vertex for large neutrino mass, Aj ^> 1 



We now consider the vertex function in the limit of large Dirac-neutrino masses (Aj 
Xj). The leading terms of the one-loop functions for Xq = Xz are:0 



Co 
C n 

C*23 
C*24 

B 1 

Co 
C u 
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C24 
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In A, 



A; 
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C 



1 In A, 



12 
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1 1 x 3 . „ , . In A 
InAi + - + -6 + Xz 



2(£>-4) 4 
+ [-30 + 5X Z + 24(4 - X z )a(y)} 



72A,. 



12Ai 
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with 



y = y/l/Xz - 1/4, 
a {y) — V arctan(l/2j/). 

The large Dirac-neutrino mass limit of the vertex function is: 



1 3 



^+ 3-^(1 



In A,; 



8c : 



IV 



+ ^ (- 66 - A ^ + 96 4^ + 5s 2 w X z ) 
+ i (-8 + 2A Z - 32 S ^ - As 2 w X z ) a{y) + /} _ , 



Numerically, this means for the value Xz = 1.286: 

1 



+ 



In A; 



V(X N ) - V(0) 



(D.31) 
(D.32) 
(D.33) 
(D.34) 
(D.35) 
(D.36) 
(D.37) 
(D.38) 

(D.39) 



(D.40) 
(D.41) 



;D.42) 



[A 4 + 2.881nA 4 - 4.47- (2.52 + 2.11 x i)] . (D.43) 



8 The functions Cn, C12, C23 do not contribute to the large mass limit of the vertex, and we 
reproduce only their first leading terms, of order In A/A. If these functions arc of relevance for an 
application, one should determine also the terms of order 1/A. 
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Here, we subtracted from the large mass limit of the vertex function its value at zero 
mass (p.66|) in order to obtain (|3.8|) . 



D.l.l A relation to the Zbb vertex 



The loop contributions to the vertex function V(X) are gauge- invariant. We calculate 
them in the 't Hooft-Feynman gauge. They describe not only the flavour- changing 
Zff vertex, but also the mass dependent terms of the flavour-diagonal vertex. There 
is one case where the effect is quite visible, namely the Zbb vertex with virtual t quark 
exchanges. The exact one-loop expression for the t quark mass dependent part W of 
the Zbb vertex correction was first given in Eqn. (22) of |37]] (form factor 5k, calculated 



in the unitary gauge). The large t quark mass limit is (Eqn. (2.4.30) of [|88 



M(Z^bb) ~ [7 M w 6 - 7^7506 + 7^(1 - 75) W(A t , A z )] tt, 
W(X t ,X z ) a 1 



a b 



V t 



nAs 2 w 



\Ptb\Vt, 



A; 



A i +(g+-^)lnA i + (9(l) 



This agrees with the leading terms in ( p.42|) 



(D.44) 
(D.45) 

(D.46) 



D.2 The vertex for small neutrino mass, A^ <C 1 

The limits of the C functions for A, = Xj and Xq = Xz 3> Aj are: 

'l + Xz + lnX z -tn)X l + ... , (D.47) 
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(D.48) 



(D.49) 



(D.50) 
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B 1 
Co 

A 1^23 
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1 1 X 
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AzC a2 = -(c -S + l)(A i -l) + ... 
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(D.51) 

(D.52) 
(D.53) 



-2(co - 5 + 1) + - [A z c - 4(c - B + 1)] A* + . . . , (D.54) 
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3-A z 

for Az = 1.286, y ~ 0.73, and the values of —Co and — Co at A, = 

Azc = y-Li2 (1^7) -^ln 2 (l + A z ) + 7rln(l + A z ) xi, 



(D.55) 
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Only the functions C develop imaginary parts, and only for 

X z > 4Ai. 



(D.58) 



(D.59) 



The functions Co,Cn,Co contain terms of the order AzhiAz at Aj <C Az, but these 
terms cancel in the form factor V. They survive in the case Az — > at constant A». 

A further note: the Euler dilogarithm is badly converging on the unit circle. In 
fact, for one of the Li 2 above it occurs |1 — Az/2 ± i x Xzy\ = 1. In that case, one may 
use: 



3?eLi 2 (e^) = 
taking the value of <ft fulfilling the condition: 

(4> ± 7r) G ( — 7T, 7r). 



^Li 2 (l) + ^(0±vr) 2 , 



;d.6o) 



(D.ei; 



The resulting small mass limit of the vertex is given in [^, |J for the case of 
sequential Dirac neutrinos: 



V(Xi^l,X z ) 



TIL ° C W 

h D-A 



+ a + a L (Xi In Xi) + a x Xi + 



(D.62) 



29 



The divergent constant was left out in the introduction, Eqn. ( |1.16[ ). 
We want to stress that in this limit, 



a L = 0. 



(D.63) 



The ansatz in Eqn. (1) of [12| allowing for 7^ is too general for the case of small 
neutrino masses in this respect. Such a term occurs for Xq = 0, see section p. 3 . 



A numerical value of ao was given in p9 |. The analytical expression is: 
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The general expression for a\ is 
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(D.65) 



with 1/ and B from ( |D.40| ) and ( p.56| ) . Further, we used Co and Co of (p.57|) and ( p.58| ) . 
With the inputs M w = 80.410 GeV and M z = 91.187 GeV, these formulae yield: 



a = 1.2584 + 1.0524 x i, 
ai = 2.5623 - 2.2950 x i. 



(D.66) 
(D.67) 



D.3 The vertex for A Q = 



The first calculation of the non-diagonal Zfif 2 vertex seems to be in ||89|| , where the 
approach was simplified considerably by the approximation Q 2 = 0. We should mention 
that this limit is of physical relevance only in the large neutrino mass limit, since only 
then it is Ai ^> Az ~ Aq = 0. If we are interested in applications where Xq > 4Aj, 
the amplitude gets essentially complex valued. The limit of f89| , however, implies 
automatically the relation Aj > Aq, even for Aj <C 1, and the amplitude is essentially 
real. So one cannot expect a continuous behaviour for the small or medium neutrino 
mass range. The virtual fermions are considered Dirac particles (Aj = Xj). 

The vertex may be easily derived from the formulae given in fllO||: 



V(X l ,X Q = 0) 
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with 



o 



'l-AJx + Ai (1-AJ 2 1-A, 



The divergent part is independent of A, and of Aq. 
Explicitly: 

_ , , . , , / A^ In A ? ' A? In A,' A? A,' 8cw 

y(A i; A Q = o) = r 3 £ 3 ./ ' +2 ' , 2 -r-V + 6 ^r-V + - 



(l-Ai) 2 (1-A,) 2 1-A, ' 1-A; D-4 

(D.71) 

with being the weak isospin of the virtual neutrinos. The constant finite terms of 
the vertex vanish for small A,: 

K(A 4 «1,A Q = 0) = J^Vn^ + e^ + ^y . (D.72) 

This has to be compared to (p.62| ) and ( p.63| ). 
The large mass limit is: 

7(A 4 »l,A o = 0) = 4 L (^ + 31nA i -5 + -^L) . (D.73) 

This has to be compared to (p.42|) and ( p.43| ). 

Finally, for the sake of completeness, the value of the vertex at the weak scale: 

y(A t = l,A Q = 0) = iJ*Q + J^Y (D.74) 
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